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Abstract 

In the first part of the paper, we introduce the concept of observable quantities 
associated with a macroinstrument measuring the density and temperature and with a 
microinstrument determining the radius of a molecule and its free path length, and also 
the relationship between these observable quantities. The concept of the number of 
degrees of freedom, which relates the observable quantities listed above, is generalized 
to the case of low temperatures. An analogy between the creation and annihilation 
operators for pairs (dimers) and the creation and annihilation operators for particles 
(molecules) is carried out. A generalization of the concept of a Bose condensate is 
introduced for classical molecules as an analog of an ideal liquid (without attraction). 
The negative pressure in the liquid is treated as holes (of exciton type) in the density 
of the Bose condensate. The phase transition gas-liquid is calculated for an ideal gas 
(without attraction). A comparison with experimental data is carried out. 

In the other part of the paper, we introduce the concept of new observable quantity, 
namely, of a pair (a dimer), as a result of attraction between the nearest neighbors. 
We treat in a new way the concepts of Boyle temperature Tb (as the temperature 
above which the dimers disappear) and of the critical temperature Tc (below which 
the trimers and clusters are formed). The equation for the Zeno line is interpreted 
as the relation describing the dependence of the temperature on the density at which 
the dimers disappear. We calculate the maximal density of the liquid and also the 
maximal density of the holes. The law of corresponding states is derived as a result 
of an observation by a macrodevice which cannot distinguish between molecules of 
distinct gases, and a comparison of theoretical and experimental data is carried out. 

1 Introduction 



When introducing the concept of observable quantity in equilibrium thermodynamics, one 
must keep in mind the fact that the observation itself should be carried out in discrete 
intervals of time that are widely separated from one another. When standing on a purely 
mathematical point of vie one must agree that the processes of establishing an equilibrium 

^In due time when the author was constructing asymptotic expansions of the Schrodinger equations in 
powers of a small parameter h, one of the presently most noted physicists told him that the asymptotics near 
the turning points cannot be considered as semiclassics because the Landau criterion for being semiclassical 
is violated there. The author, as a mathematician, believed that the asymptotics even in the domain of deep 
shadow and the "instantons" obtained for the imaginary number h can still be considered as the semiclassical 
asymptotics. Recently, Yu. M. Kagan clearly explained the author that the physicists mean only the case 
fj. < when speaking about the Bose-Einstein distribution. But the author, as a mathematician, believed 



1 



require infinite time. However, in mathematics there are some concepts which are similar to 
the notion of "half-life" in physics. For example, one can introduce a time interval during 
which the difference between the current state and the state of equilibrium in the course of 
relaxation becomes e times less. 

In approximation theory and in the theory of numerical methods, especially after the 
well-known paper of Mandel'shtam and Leontovich [39] , the following relaxation process was 
in use: at first, a reacting system is brought to some equilibrium. Then one rapidly changes 
one of the conditions (e.g., the temperature or the pressure) and traces the evolution of the 
system towards a new equilibrium (see, for example, the article "method - relaxation" in 
the Great Encyclopedia of Oil and Gas, http://www.ngpedia.ru [in Russian]). 

Since the observation intervals should be "equal" to the relaxation time, they are large 
enough, and one can refer to the process as the multi-step relaxation process (MRP). Eco- 
nomic and historical processes, and also biological processes in a living organism, belong 
to phenomena of this kind, and therefore, from time to time, thermodynamic models of 
these processes arise. The formation of clusters, according to the scheme suggested below in 
Sec. 3.2, can serve as an example of a multi-step relaxation process. 

The fact that time intervals of observation are discrete is the most important point to be 
taken into account when speaking about the instruments of observation. 

The difference between readings of measuring macro- and microinstruments in thermo- 
dynamics is related to the following aspects. 

1. A macroinstrument does not take into account the motions of nuclei, of electrons, 
and even of atoms within a molecule and regards any molecule as an individual particle. 
Mathematically, this corresponds to imposing rigid constraints on the elements forming the 
molecule. That is, we must modify those axioms of mechanics in which we consider all 
elementary particles and their behavior in the configuration space whose dimension is equal 
to the tripled number of elementary particles. 

2. A macroinstrument measuring density counts the number of particles in a fragment 
of the volume; however, it cannot trace the movements of particles with different numbers 
during discrete finite time intervals. At each discrete time moment, this device counts the 
number of particles in the same fragment; however, it cannot notice what is the exact position 
of any particle indexed at the previous time moment and whether or not this particle really 
is within the chosen fragment. Mathematically, this means that the arithmetical law of 
rearrangement of summands holds. The sum does not depend on the way in which we have 
indexed the particles. In this sense, the laws of classical mechanics are even modified in a 
more substantial way. 

Let us quote from the textbook [1] on quantum mechanics, where the authors define 
the basic property of classical mechanics: "In classical mechanics, identical particles (e.g., 
electrons) do not lose their 'personality,' despite the identity of their physical properties. 
Specifically, you can imagine that the particles forming a given physical system are 'indexed' 
at some time moment and then one can trace the motion of each of the particles along its 
own trajectory; then the particles can be identified at any time moment. ... In quantum 
mechanics, it is fundamental that there is no way to trace each of the identical particles 
separately and thus to distinguish them. We can say that, in quantum mechanics, identical 

that this restriction is artificial and considered system (l)-(4) in the general case, without any restrictions 
on the number iV^ of particles at the ith energy level. But the natural restriction J2 — ^ ' — ^^ill 
exists and is taken into account by the author. The Bose-Einstein condensate also exists but in a small 
neighborhood of the zcroth energy level (small compared with n) rather than at a single point. The author 
continues to use the name "Bose-Einstein" for the obtained distribution and the condensate phenomenon. 
The general asymptotics is constructed for N :$> log A^, and this asymptotics holds for N = 100. 
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particles completely lose their individuality" (Russian p. 252). 

A macroinstrument does not keep this basic property either. Mathematically, this means 
that, to take this property into account, we should impose some new constraints, which are 
already explicitly nonholonomic, on the mechanics of many particles and, which is especially 
important, we should take into account the permutability of particles in the definition of 
density, namely, any permutation of particles does not modify the density. 

In thermodynamics, the gas molecule density is measured. Although the gas molecules 
differ from each other and the Boltzman approach to studying the molecules is consistent 
with the objective reality, the difference between the molecules does not play any role when 
the molecule density is determined. If the density is considered in a small fragment of the 
vessel, which contains approximately a million of particles, then it turns out that the density 
in this fragment coincides with the average density in the entire vessel up to a thousand of 
particles (up to 0.1 %) and is independent of the particle numeration. 

It follows from these considerations that the entropy (in contrast to the Boltzmann- 
Shannon entropy) should take into account the permutability of the indices of the particles 
(cf. [2], Sec. 40). 

Hence, for an ideal gas 

j j j 3 

-^{S + aN + I3E) = (2) 

any- 
where Hj stands for the average number of particles in each of the Gi states of the jth group 
and a and (3 are some constants (see [2], the footnote on p. 184, and also jl] and [5]), the 
entropy must be of the form 

S = Y.^{G, + N,) log(G, + N,) - N, log N, - G, log G,}, (3) 

3 

S = ^G,[{1 + nj)\og{l + n^) - log n,] . (4) 

3 

In other words, the entropy has exactly the same form as in the Bose-Einstein quantum 
case. This face is proved for balls and boxes in in the footnote in Sec. 46; also see [H [13]. 

We have noted above that a macroinstrument and its measurements force mathematicians 
to reorganize even the axioms of classical mechanics. However, mathematicians are forced 
to do so by entering the corresponding small parameters and passing to the related limits. A 
macroinstrument and its measurements still reduce the time spent to perform constructions 
of this kind. However, when one speaks of the axioms of thermodynamics, which is based 
on laws derived by great physicists who used ancient experiments conducted on Earth, it 
then turns out that the above considerations modify the classical concept of thermodynamics 
completely. Meanwhile, microinstrument^ also play a role in classical thermodynamics; they 
enable one to calculate the dimension related to the number of atoms in the molecule. 

In the mathematical literature, as a rule, the number of degrees of freedom coincides with 
the number of independent generalized coordinates. However, there notions are distinct in 

^In mathematics and mechanics, the difference between micro- and macro-observations is defined as 
foUows: "the radius a of a molecule is much less than the typical length of the vessel (provided that the 
shape of the vessel is given)", i.e., there are two scales in the problem, which correspond to macro- and 
microinstruments . 
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the standard thermodynamics, because the volume is three-dimensional, which is established 
by the macrodevice, whereas the number of degrees of freedom is related to the number of 
atoms in a molecule and is measured by the microdevice. 

Let us explain the following experimental fact. In some cases, the number of degrees 
of freedom for diatomic and polyatomic molecules is an integer. In our opinion, this hap- 
pens because the intramolecular communications (the distances between the atoms of the 
molecule) are very hard, and, when the temperature increases, no new degrees of freedom 
arise. Generally speaking, the number of degrees of freedom fundamentally depends on the 
energy of the molecules, and the energy of different molecules of the same gas is different, 
and, apparently, to the average energy (the temperature) there must correspond the average 
number of degrees of freedom, which is hence must be noninteger. However, on one hand, 
tight connections enable one to excite almost all molecules for a sufficiently high (room) tem- 
perature and, on the other hand, to give the molecules no possibility to excite new degrees 
of freedom (e.g., the vibrational ones). If the connections are not so rigid, then the number 
of degrees of freedom depends on temperature, and hence on energy, and is not an integer 
in general. This is clear from the comparison of the values of the heat capacity Cy with the 
experiment: for hydrogen sulfide with three atoms, the theory gives 5.96, and the experi- 
ment 6.08, for carbon dioxide, the experiment gives a greater value Cy = 6.75 (T = 15C, 
P = latm), and, for carbon disulfide, the vale is almost two times larger, namely, 9.77. In 
the case of diatomic molecules, say, for nitrogen, the theory gives 4.967 and the experiment 
shows 4.93; for the chlorine, the value is almost 20% higher, namely, 5.93, etc. 

It turns out that the number of degrees of freedom coincides with the dimension of the 
generalized Bose gas which is regarded as a distribution of a classical gas. 

Landau and Lifshitz notice this fact for the three-dimensional Bose gas. They write 
that these equations (PV^^^ = const) coincide with the equations of the adiabatic line for 
an ordinary monatomic gas. "However, we stress," the authors write further, "that the 
exponents in the formulas VT^^"^ = const and PV^^^ = const are not related now to the 
ratio of specific heat capacities (since the relations Cp/cy = 5/3 and Cp — = 1 fail to 
hold)" [2], p. 187. 

One can show in a quite similar way that, for the five- dimensional and six- dimensional 
Bose gas, the "Poisson adiabatic line" coincides with the Poisson adiabat for the two-atomic 
and three-atomic molecule (see [2], Sec. 47, Diatomic gas with molecules of different atoms. 
Rotation of molecules). With regard to the above stipulation, as — )■ — oo, we obtain 
precisely both the condition Cp — = 1 and the ratio Cp/cy coinciding with relations well 
known in the old thermodynamics. 

Remark 1 . The three-dimensional case of the Bose-Einstein-type distribution can be repre- 
sented as 

i+k+m=j 

i,k,m j i+k-\-m=j 

Y j^i,k,m = Y^ Y ^i,k,m = Y^^r 

j i+k+m=j j i+k+m=j j 

The Bose-Einstein "average" values nj^fc^m of the occupation numbers A^j,A;,m depend only on 
the energy, i.e., on the sum i + k + and 
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so that 

- Qj (j + 2)! 

3 nt,k,m g/3(i+A;+m-/.) _ X ' j!3! ' 

The transition to integer dimensions is similar; the fractional dimensions are obtained by 
passing from factorials to F-functions. A more rigorous approach in described in [T3l H2l H3] . 

2 A new ideal gas and a new ideal liquid 
as observable quantities 

2.1 The number of degrees of freedom for T <Tc and P < Pc 

Let us now proceed with finding the number of degrees of freedom for for a low temperature 
that does not exceed the critical one: T < T^. 

The Maxwell-Boltzmann equation for the ideal gas is of the form 

PV = NT, (5) 

where P stands for the pressure, V for the volume, N for the number of particles, and T for 
the temperature. 

Denote by Z the dimensionless quantity Z = j^, which is called the compressibility 
factor. Equation (5) can be represented in the form Z = 1. Let us express the Bose- 
Einstein-type distribution for the fractional dimension D using polylogarithms. 

Represent the thermodynamic potential of the Rose gas of the fractional dimension D = 
2(1 + 7) the form 

7^2+7 poo y.1+7 JJ. 

n(/z,T) = (Cm)i+^r^;^^y^ j-^j^ = -T'+\CmY^^'VU,^,{a), (6) 

where T stands for the temperature, m for the mass, C is a constant, a = exp{fi/T) is the 
activity, /i is the chemical potential, and F stands for the Euler gamma function . 

The function Lis(a) introduced in (6) is referred to as a polylogarithm and is defined by 
the rule 

-| i>oo 

where ({s) stands for the Riemann zeta function. 

To pass to the dimensionless units, we introduce the temperature Tr in such a way that 
T = TrT,. 

The expressions for the dimensionless pressure Pr = P/ Pc and for the number of particles 
N that correspond to the thermodynamic potential (6) are of the form 

T^^Li,+,(a) 
= C(2 + 7) ' 

iV = \/T;+^Lii+^(a). (9) 
We have (for the definition of 7c, see below) 

Q' = -T^+^{Cmy^-WU2+^{a). (10) 



The following formula can thus be obtained for the compressibility factor Z: 



^ ^ Li2+^(a) 
Lii+7(a) ' 

In particular, for a = 1 (i.e., for yU = 0), we have 

07 + 2) 
C(7+l)' 



fill 



(12) 



As is well known, in the Bose-Einstein theory, the value /i = corresponds to the so- 
called degeneration of the Bose gas. 

For a classical gas satisfying the same relations, the degeneration coincides with the 
critical point T = T^, P = Pc, and Z = Z^.- Consequently, one can write 7 = 7c for Z = Z^. 
in (12), namely, 

= ^!2£±4. (13, 

and to every pure classical gas there corresponds its own value of 7c. 

The entropy in the dimension D = 27 + 2 can be evaluated in the standard way. The 
great thermodynamical potential is considered, 

VT 1 t^+^ dt _ -VT^+^ 



A2(i+7) r(2 + 7) Jo (eVa) - 1 A'2(i+7) 

where A' = const/ (271771)^/^, the dimension D is equal to 27+2, T stands for the temperature, 
and a = e^^^ for the activity. 
The number of particles is 



The compressibility factor is 



Let us evaluate the entropy. 



PV _ Li2+7(a) 



NT Lii+^(a) 



(2 + 7) Lii+^(a) - Lii+^(a)^ 



A'2(l+7) 

For /i = 0, = T/Tc, and Pr = P/ Pc, the specific entropy is equal to 

|Ut.^i = (2 + 7)C(2 + 7). (14) 

E. M. Apfel'baum [Apfelbaum] and V. S. Vorob'ev ^ compared the Bose distributions 
of fractional dimension D = 27c + 2 in the (P, V) diagram with the experimental critical 
isotherms for various gases. We present these graphs in Figs. 1-5. In Fig. 5A, the graphs for 
the nitrogen and oxygen are shown, which have been constructed by Professor V. S. Vorob'ev. 
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Figure 1: (a) Isotherms of pressure for the van der Waals equation are shown by sohd hues. 
The hues formed by circles are constructed from computations for 7 = 0.312 (i.e., for the 
ideal "Bose gas"), Z„ = 3/8, p = P/Pc, and n = N/N^. 

(b) Isobars of density for the van der Waals equation are shown by solid lines. Line 1 is the 
binodal. The circles correspond to isobars of the "Bose gas" for 7 = 0.312. 




Figure 2: Critical isotherms for the Lennard- Jones system. Symbols 1 and 2 correspond to 
numerical calculations. Line 3 corresponds to the ideal Bose gas for 7 = 0.24. 
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Figure 3: (a) Isotherms for argon. The sohd hnes correspond to experimental data. The hne 
formed by circles is constructed in accordance with the isotherm of the ideal "Bose gas"; 
Zcr = fgg} = 0.29, p = P/P,, and n = N/N,. 

(b) The same for water, Zcr = 0.23. 

(c) The same for copper, Zcr = 0.39. 




Figure 4: Isotherms for water. Symbols 1 and 2 correspond to experimental data, and line 3 
corresponds to the computation for the Bose gas. 
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Figure 5: Critical isotherms for mercury. Symbols 1 and 2 correspond to experimental data, 
and line 3 corresponds to the computation for the Bose gas. 




Figure 5 A: 1. Isotherms for nitrogen. 2. Isotherms for oxygen. 



2.2 Bose condensate as an observable quantity 

in classical thermodynamics. Relativity principle for MRP 

We shall show that the Bose condensate in classical thermodynamics is the condensate of gas 
(vapor) into liquid (in contrast to the statement presented in the manual |2j in the footnote 
on p. 199). 

Example 1 . Consider the example given by a famous theorem in number theory, namely, 
the solution of an ancient problem, which has the Latin title "partitio numerorum." This 
task involves an integer M which is decomposed into N terms, for example, if M = 5 and 
N = 2, then 5 = 1 + 4 = 2 + 3, which gives two solutions to the problem, Ai = 2. 

If M = 10^^ and = 1, then the decomposition has only one version, and Ai = 1. If 
M = 10^^ and = 10^'^, then there is also only one version of decomposition, namely, the 
sum of ones, i.e., Ai = 1. 

Obviously, there is a number Nf. for a fixed M such that the number of versions of the 
decomposition, Ai, is maximal possible (this number is not unique in general). The number 
log2 Ai is referred to as the Hartley entropy. At the point at which Ai reaches its maximum, 
there is a maximum entropy. 

Let a partition M = ai + ■ ■ ■ + of M into summands be given. Denote by Nj the 
number of summands on the right-hand side that are precisely equal to j. 

Then the total number of summands is Ylj this number is equal to A^, since we 

know that the total number of summands is A^. Further, the sum of the parts equal to j 
is jNj, since there are Nj summands, and then the sum of all summands can be obtained 
by summing these expressions over all possible j, i.e., Yljj-^j^ ^^'^ ^^i^ is equal to M. 
Namely, 

oo oo 

J2n, = N, Y.'N, = M. (15) 

i=l i=l 

The very nonuniqueness of the above maximum and an uncertainty concerning the num- 
ber of the maxima enabled Erdos to obtain results with accuracy up to o(a/M) only. 
Thus, the Erdos theorem holds for the system of two Diophantine equations 

oo oo 

Y^N, = N, Y.'Ni = M. (16) 

i=l i=l 

The maximum number of solutions of the system is achieved provided that the following 
relation holds: 

AT, = /3-^M]/2 log Me + aMl'^ + o{Ml''^), (3 = tt^/2/3, (17) 

and the coefficient a is defined by the formula (3/2 = e""^/^. 

If one increases the number A^ in problem (16) and keeps the number M constant, then 
the number of solutions decreases. If the sums in (16) are counted from zero rather than 
from one, i.e., 

oo oo 

J2^Ni = {M-N), ^A^i = Ar, (18) 

i=0 i=0 

then the number of solutions does not decrease and remains constant. 

Let us explain this fact. The Erdos-Lehner problem [6] is to decompose a number 
into N < Nc summands. 

The decomposition of the number 5 into two summands has two versions. If we include 
also 0, then we obtain three versions, 5+0 = 3+2 = 4+1. Thus, the inclusion of zero gives 
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the opportunity to say that we decompose the number into k < n summands. Indeed, the 
expansion of the number 5 into 3 summands includes all previous versions: 5 +0 +0, 3 +2 
+0. and 4+1 +0 and adds new options that do not contain zero. 

Here the maximum does not change much [6]; however, the number of options cannot 
decrease, because the zeros enable the maximum to remain constant, and the entropy never 
decreases; after reaching the maximum, it becomes constant. This very remarkable property 
of entropy enables us to construct a general unbounded probability theory [7]. In physics, 
the effect is identical to the so-called phenomenon of Bose condensate. 

Let us pose the following question: what is the difference between arithmetic, together 
with the problem of "partitio numerorum," and the Boltzmann-Shannon statistics? If we 
assume that 4 + 1 and 1 + 4 are two different versions, then we obtain the Boltzmann- 
Shannon statistics. The number of versions of decomposition, Ai, is growing rapidly. Thus, 
the "noncommutativity" of the addition gives additionally a huge number of versions of 
decomposition, and the Hartley entropy (which is equal to the logarithm of the number of 
versions) coincides with the Boltzmann-Shannon entropy. 

Therefore, we have proved that, if we add zero to the family of possible summands and 
decompose a number M into summands, then this is equivalent to solving equations (18), 
i.e., to imposing relations for the number of particles and for energy. Here the number of 
zeros increases drastically; if M = 5, then, for > 3, the number of zeros is 22. However, 
the number of ones is also large, although it is twice smaller than the number of zeros. 

It is very visible to consider the Bose condensate as the number of zeros; however, this is 
inaccurate. The Bose condensate occurs in a neighborhood of a point at which the energy 
vanishes rather than at the point itself. Nevertheless, if one writes poS{k.) (where po stands 
for the density and the vector k is the momentum) for the Bose condensate at rest, then 
this notation is true, because the density is the limit 

Nm 

Po = hm —— , 

where stands for the number of particles, V for the volume, and m for the mass of the 
particle. This means that, as — ?■ oo, the bell-shaped function near the zero energy is 
converted to the 6 function. 

By an ideal (or perfect) liquid we mean a liquid without attraction and without any 
surface tension. This is a liquid which can exist for a positive pressure only together with a 
saturated steam. 

In this case, the perfect liquid is the result of an optical illusion, and this "perfect liquid" 
is the same ideal gas in the condensate, with another density. It can be described, as in the 
case of a consideration of a Bose condensate, in the form po^(k); where po stands for the 
density of the condensate. It cannot exist without a volume trap, which is similar to the 
case in which a container with gas has a hole, and there is a vacuum outside the vessel, in 
which case this liquid, which looks as if it is boiled, is going away together with the gas. 
The mean speed of the particles inside the liquid is the same as the mean speed of the gas 
particles. This corresponds to the condition that the temperature in the system "liquid- 
saturated steam" is the same in the liquid and in the gas. The liquid in a closed vessel is 
a fluctuation standing at a fixed place (5(k)), or, speaking in a simpler way, this liquid is a 
"resting" Bose condensate (cf. [8], p. 204). 

Small crystals that occur in a supersaturated solution coincide with the Bose condensate 
only if they are not composed of mutually connected particles, and, moreover, if the particles 
are continuously exchanged with the particles in the solution; moreover, the small crystals. 



11 



as solids, are an optical illusion, namely, we simply do not see that the particles of crystals 
are permanently transposed with particles of the solution. In other words, this is by no 
means a crystal, this is a fluctuation; however, this fluctuation is relatively immobile. 

Thus, the Bose condensate for classical particles represents some "special density fluctu- 
ations;" only this, and nothing more. 

A. I. Anselm constructed his theory starting from the Eyring formula for free energy. 
The liquid structure model accepted by Eyring is in fact closer to the strongly compresses 
gas model [H, 05]. 

One can talk about the density in a "special cluster fluctuation" of a part of our vessel 
with gas. If we speak of the density in this cluster only, this means that (as in the example of 
a small volume with one million of particles) one cannot speak of the number of particles that 
are placed in the cluster as if they are frozen and do not move. This is only an appearance, 
and all of the particles or a part of them can be replaced in a minute by another ones, and 
the indexing inside the "fluctuation" cluster can change every minute. At the next time 
step, this can be the same picture but probably with different particles involved. 

We speak about some fragment of the volume. In fact, the particles that are more 
concentrated can be spread out over the entire vessel. However, if there is at least a little 
gravity of the Earth, then the fluctuations with more concentrated particles accumulate 
near the bottom. If we consider a vessel with gas in the form of a perfectly reflecting sphere 
(see [9]-[12]), then, due to the repulsive force occurring at the border, fluctuations of this 
kind are located near the center of the ball. 

From the standpoint of our observation, in discrete time intervals at far distances from 
each other, the denser fragment of the volume, i.e., the Bose condensate, is at rest and hence 
corresponds to a small momentum in the Bose-Einstein-type distribution. Mathematically, 
the MRP model corresponds to this phenomenon. This property will be called the relativity 
principle for MRP. 

Let us repeat once again that the only fact which can be guaranteed by the generalized 
theory of Bose condensate is that there will be a higher density of particles at the bottom. 

Example 2. Let a gas be contained in a closed vessel at a room temperature, and let 
the gas be almost satisfying the Clausius relation 

P = pT. (19) 

We cool the vessel down to a temperature T = 0. At some temperature T > To, a liquid 
is formed. The temperature Tq is referred to as the dew point. According to the standard 
conception, the fluctuations above the temperature of the dew point are of the order of a/ZV. 
After the formation of liquid, the gas, which is called a saturated steam in the physical 
literature, also satisfies relation (19). It is quite rarefied. According to the van der Waals 
model, there are no singularities at the dew point under the gas-liquid passage (on the so- 
called binodal). According to experimental data, there are no large fluctuations either in the 
usual sense at the dew point. 

Finally, the most important thing. The experiment shows that, at T = Tq, the gas is 
rarefied, and it remains an ideal gas in the sense of relation (19), i.e., in the Boltzmann- 
Maxwell sense. 

There is, however, a fluctuation of the type of a stationary Bose-Einstein condensate, in 
this fluctuation, the molecules by themselves placed inside this fluctuational fragment can 
possibly move with the same velocities as those of the gas molecules and, if it were possible 
to enumerate them, then the numbers will be changed very quickly. If shall refer to this 
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fluctuation (of tlie form of the Bose-Einstein condensation) as liquid, then actual molecules 
of the liquid move in it with the same speeds as the gas molecules (of the "saturated vapor"). 

To represent this picture in a more visible way, imagine a bunting which winds from one 
roller to the other. Between the rollers, under the material, a strong wind blows from a hose. 
We see a " hump" is formed between the rollers; however, it can be assumed that we do not 
see that the bunting moves. 

Nevertheless, as the density of the Bose condensate increases, our macroinstrument can 
fix the bound of the density and show us that there is a more dense phase and a less dense 
phase. Hence, only the original macroinstrument can show us the bound of this abstract 
liquid, i.e., of the second phase. 

First, the Bose-condensate at rest, i.e., the gas compaction, is being formed (because of 
the relativity principle for MRP), and then there arise quantum forces, i.e., attraction forces 
(see below), acting on the "nearest neighborhoods", the more so because the molecules move 
slower at a low temperature. 

If liquid droplets occur below the temperature of the "dew point," then the droplets 
are spherical, even under the presence of the gravity of the Earth (physicists refer to the 
very gravity, as a rule, when claiming that the border between gas and liquid is fiat). The 
pressures in the droplet and in the gas (the saturated vapor) are different, due to the surface 
tension. 

Therefore, the main rule of the equilibrium "vapor-liquid," namely, the coincidence of 
the of pressures, really holds at the dew point only if we neglect the surface tension, and thus 
neglect the attraction of liquid molecules, because these two effects are inseparably linked 
with each other. Our concept of a new ideal gas is based on the very assumption on the 
absence of attraction between the molecules. 

The picture in which the attraction and the surface tension play no role can be is observed 
in experiments if the temperature is equal to the temperature Tq of the gas-liquid transition 
(i.e., at a point of the "binodal") and Tq is still greater than the temperature at which a 
droplet of critical radius has been already formed. Then, at T = Tq, the incipient drops 
spontaneously shrink and occur at another point. These drops cannot live without the 
surrounding saturated vapor; one can see these drops but cannot feel them. 

If the vessel is spherical and the mean free path is comparable to the size of the vessel 
(similar to the so-called Knudsen criterion; see [9]-[T2], then the probability of such a 
virtual drop is larger at the center of the vessel. 

In this case, if we make the labelling of several molecules by launching few isotopes which 
can be traced, then these isotopes will pass freely through the liquid to vapor and back, and 
they will form a a denser structure near the center of the ball, in such a way that, when 
illuminated by parallel rays, it will provide a shade. However, it is impossible to take this 
drop from the gas medium. One can see an ideal liquid but cannot feel ilH. Possibly it is 
better to refer to it "virtual liquid." 

This approach is unusual for the majority of physicists. Although everyone knows that, 
say, when photons are collected at a focus at which their "density" is high, then it is impos- 
sible to separate the focus from the "photon medium." 

A mathematical analog of the quantum Bose condensate for a classical gas is a liquid 
without attraction in which the speeds of the molecules are approximately the same as the 
speeds of the molecules in the saturated vapor. The attraction between molecules results 
in a significant correction provided that the radius of a drop is greater than the critical 
one; however, this correction abolishes the conditions of the vapor-liquid equilibrium for the 

''One cannot drink it but can breathe it in. 



pressure. Therefore, the problem must be divided into two separate problems, namely, 1) an 
ideal gas and a perfect liquid without attraction, and 2) the consideration of the attraction 
for the case in which the decay into two phases has already been carried out and the radius 
of the drop exceeds the critical value. 

Remark 2. We define the temperature from the overcondensate part of the system, i.e., from 
the gas until the volume of particles in the condensate is comparatively small, i.e., until the 
surface tension is formed and a drop of critical radius size appears. 

The drop of critical radius size is the result of a different MR- process, i.e., of the quan- 
tum dipole-dipole interaction clusterization according to the scheme given in Sec. 3.2. The 
nucleation process consists of two mutually related MR-processes. The Bose condensate in 
the first MR-process is the nucleation starting mechanism including the quantum effect of 
dipole-dipole attraction and the quantum effect of exchange interaction of identical particles. 



2.3 Asymptotic continuation of a perfect liquid 

to the second sheet as the volume of the liquid increases 

In the manual by Landau and Lifshitz and in other manuals, the spectrum is calculated by 
the Weyl-Courant formula. Such calculations require the use of the phase volume, and the 
volume V of the configuration space naturally arises. We determine the spectrum starting 
from the number of degrees of freedom and actually use the volume only in the final result 
to pass from the number of particles to the density. As was already seen, the number of 
degrees is equal to the dimension of the Bose-Einstein-type distribution. 

The gas spinodal, which is defined in a new way as the locus of isotherms of a new ideal 
gas, is formed at the maximum entropy at the points at which the chemical potential /i 
vanishes. 

Therefore, on the diagram (Z, P^), the spinodal is a segment < 1, Z = in the case 
of the van der Waals normahzation Tj. — T/T^ and Pr = P/Pc- 

Until now we, maximally following the traditional notation used in [3]. preserve the 
volume V, although neither the equation for the ^.-potential given in [3, $ 28] 

dQ = -SdT-Ndii (20) 

nor relations (l)-(4) contain the volume V . We interpret the Bose-Einstein condensate as 
a liquid phase, and because for N > the number of overcondensate particles remains 
constant, the hquid is "incompressible". 

For Tr < 1, the Bose condensate occurs and, consequently, for the liquid phase on the 
spinodal, the quantity 

iV = T>+iC(7c+l) 

remains constant on the liquid isotherm. This means that the isotherm of the liquid phase 
that corresponds to a temperature is given by 

P P 

- ' (21) 



TrN r;^+^c(7c + i) 

All isotherms of the liquid phase (including the critical isotherm at = 1) pass through 
the origin Z = 0, P,. = and then fall into the negative region (or to the second sheet). The 
point Z = corresponds to the parameter 7 = 0, and hence to the continuation to 7 < 0, 
since, for = 0, the pressure 
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can be continued to > 7 > — 1. 

We shall see below that the value of Z as /i — )■ is also positive, and therefore the 
spinodal for > 7 > — 1 gives the second sheet on the diagram {Z, P); it is more convenient 
to map this sheet onto the negative quadrant. 

Under the assumption that the transition to the liquid phase is not carried out for = 1, 
we equate the chemical potentials fi and /i for the "liquid" and "gaseous" phase on the 
isotherm = 1 (this fact is proved below). 

After this, we find the value of the chemical potential corresponding to the transition 
to the "liquid" phase for < 1 by equating the chemical potentials of the "liquid" and 
"gaseous" phases. 

In this section, we find the point of the isotherm-isochore of the liquid as the quantity 
X = — /i/T tends to zero. 

First of all, we take into account the fact that Nc is finite, although it is large, and hence 
we must use the obtained correction. 

In fact, the transition to integral (6) from the integral over momenta in by using 
the replacement p^/2m = e corresponds to the transition to the energy oscillatory "repre- 
sentation" or, which is the same, to the natural series. The differential de means that the 
discrete series must be taken with the same series in e, and this is precisely the natural series 
multiplied by a small parameter. 

Historically, such a representation was present already in the initial Plank distribution. 
The transition from the discrete representation of the natural series to the integral repre- 
sentation will be described in this section. This representation associates the Bose-Einstein 
distribution with the number theory considered in Example 1. On the other hand, it stresses 
that the discrete Bose-Einstein-Plank distribution depends only on the number of degrees 
of freedom and is independent of the three-dimensional volume V. 

Obviously, the discrete decompositions leading to integral (6) are not unique. Usually, 
the physicists reduce discrete decompositions to integrals over momenta and try to relate 
them to the volume V (and the phase volume, respectively). Using the natural series and 
the parameter 7, we thus stress the difference between these approaches. 

Let us construct the thermodynamics of the ideal Rose gas with boundedly many states 
at a given quantum level. Since Ni < N because of the left equality in formula (11-1), this 
condition cannot be an additional restriction. Summing the finite geometric progression, we 
obtain 



For the number of particles, we have the formula N = —dQ/dfi (see (20)). Omitting the 
volume V, we obtain 




(23) 



The potential f2 is equal to the sum i7j over i: 




(24) 



N 



A2(l+7) 



1 




exp (-^ + i) - 1 exp [{N + 1)(-^ + i)] - 1 



{N + 1)P 



(25) 



The volume V in relations (j23|) was required only for the normalization, for the transition 
form the number to the density. For 7 > 0, it does not interfere with the asymptotics as 

— )■ 00, because the term containing + 1 in the right-hand side is small. At the same 
time, it agrees with the pressure, because P = —dfl/dV. 

For 7 < 0, we omit the volume V, because even for 7 = due to Example 1, there appears 
a term of the form logA^ which must be taken into accoun10, because we have log A^ ~ 15 in 
the two-dimensional case. 

In the two-dimensional trap, the number A^ is significantly less, but even for A^ = 100, 
logA^ = 2, we can use the asymptotic formulas given below. 

On the other hand, the relation between thermodynamic parameters allows us to decrease 
the number of independent variables from three to two (cf. Fig. 9 in the variables p, T and 
Figs. 11-16 in the variables Z,P). 

Estimates. 

□ Taking the parameter 7 into account we use the Euler-Maclaurin formula to obtain 
V (-J1 4>-) = i r , J ) i.' + R, (26) 

j ° 

where = ^ + 1, k = N + l,b = 1/T, and x = —fi/T. Here the remainder R satisfies the 
estimate 

1 /"^ , „ ^ , 1 k 



1 1"^ 

R\<- \f'{x)\dx'^, where f{x 
a Jo 



^bx+x ]^ Qk{bx+>c) 

We calculate the derivative and obtain 

^ ~ (^^kibx+K) — 1)2 ~ (^gbx+x _ ]^)2' 



(27) 



We also have 



+ il){z)^ where il){z) is a smooth function and \il>{z)\ < C(l + \z\) ^. 



(e^ - 1)2 ^2 

By setting z = y and z = ky, we obtain the estimate for R: 

1 /"°° 

\R\<-n; \^{k{y + x))-i,{y + x)\dy'' 



7 -a poo -1 



with a certain constant C. For example, if x ~ (logA;)-^/^ then \R\ preserves the 
estimate|/2| ~ 0{b~"). □ 



^In this example where D = 2 and 7 = 0, there is no area ©. And this confuses speciahsts in ther- 
modynamics. Indeed, on one hand, N/G — > const, but on the other hand, it fohows from (17) that 
log Mc ~ 2 log iVc, and hence, by (17), the limit of N/& as iVc — > 00 and 6 — > 00 tends to infinity. 
This finally leads to a false conclusion that the Bose-condensate exists only for T = in the two-dimensional 
case. In fact, it exists for Td = "^^^ ( ) 15^77 ' ^^"^ ^ very small value (see Corollary 1 below). 



The energy will be now denoted by M, because without multiplication by the volume V, 
this is not the usual thermodynamics but rather a certain analog of the number theory (see 
Example 1). 

Taking account of the fact that, for the value of M, the correction in (!23l) can be neglected 
for the value of M, we obtain 

~ «r(7 + 2) j ~ 'b^ Jo i^^' ^^^^ 

where a = 7 + 1, 6 = l/T^. Therefore, 



Mi/(i+«) V«r(7 + 2) Jo e€ - 1 

We obtain (see |18]) 

y (—11 = iri^ ^ I dc + o(b~n 

Ve^J+- - 1 e^'^^+'^'^J a Jo \ - 1 e'=''« - 1 ( ^ ^ ^ ' 



ah 

= ^(A:^-"-i) + o(r"). 

By setting k = N\fi/T=o ^ 1, we finally obtain 



f°° /I 1 \ 

Ar|^/^=o = (A^^-M7))'/^'+^^7^r, where 0(7)=/ b-;^)^"^^. (30) 



Corollary 1. [Erdos formula] It can be proved that x — > gives the number N with satis- 
factory accuracy. Hence, 

Consider the value of the integral {with the same integrand) taken from e to 00 and then 
pass to the limit as e ^ 0. After making the change bx = ^ in the first term and bN^x = ^ 
in the second term, we obtain 





poo 


d^ 


b. 


Lb 






rebNc 




b. 


Jeb 





V^ + o(r^)4P^ 



+ 0{b-') = T -r^ + 0{b-') (31) 



ebNc ^ ^ '-'Jeb 

0{b~') = ^{log(£6iV,) - log(£6)} + 0(6-1) ^ 1 i^g^^ ^ Q^^-iy (32) 
b b 

On the other hand, making the change bx = ^ in ( |29|1 . we obtain 

1 r 



62 

This gives 



M. 



'M/^ll -^y, A^.= log M(l + 0(1)) + 0(6-1). (33) 
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Now let us find the next term of the asymptotics by setting 

= c-^M^/^ log c-^M^/^ ^ ^^1/2 ^ o{M^/^), where c = ^ . 

v6 

Furthermore, using the formula 

iV, = log ^ 0{b-^) 



and expanding log in 



we obtain 



a 



c-ilogc-iMi/2 



c 

a = -2 log-. 



Thus, we have obtained the Erdos formula l4^ . 

The relation N = T^''~^^({'yc + '^) is consistent with the linear relation = A{'y)Tr, where 
A(7) = (A>-^c(7))i/(i+^), for Pr < 0. 

We can normalize the activity a at the point Tc, and we can find oq by matching the 
liquid and gaseous branches at Tc for the pressure Pr = 1, in order to prevent the phase 
transition on the critical isotherm at = 1. 

In what follows, we normalize the activity for < 1 with respect to the value of ao 
computed below. Then the chemical potentials (in thermodynamics, the thermodynamic 
Gibbs potentials for the liquid and gaseous branches) coincide, and therefore there can be 
no "gas-liquid" phase transition at T,. = 1. 

Now, for the isochore-isotherm of the "incompressible liquid" to take place, we must 
construct it with regard to the relation N^. = ({•jc + 1), i.e., 

iV(T,)=T>+iC(7c+l). 
We obtain the value jiTj.) from the implicit equation 

A(7) = r>c(7c + i). 

Thus, for each < 1 we find the spinodal curve (i.e., the points at which /t = 0) in the 
domain of negative 7 |18] . 

A(7-7c)/{l+7)c(^)l/(l+7) = T7t(7c + 1), (34) 

In the set of two values of 7 corresponding to the solution f l34p . we choose the value 
associated with the largest entropy, i.e., the quantity largest in absolute value and denote it 
by 7(Tr). For = 1, we choose the value of A so that both solutions 7(1) coincide, and we 
write 7o = 7(1). 

Let Ug = e~^/^ be the activity of the gas, and let a/ = e~^/^ be the activity of the liquid. 
W e present the condition for the coincidence of M and of the activities at the point of the 
phase transition: 

= Al-(^^)l+>T-l-(^'-)l Li,_|,(^,)| (^) , (35) 
Li2+^o(ao) = 1, ag = —. (36) 



C(2 + 7c) ^"""^^^ ' ' ao' 



18 



Definition 1. The relation ag = ai/ao will be called the normalization of activity on the 
critical isotherm. 

Relations (I35l)- (j36ll determine the value of the chemical potential fi = fi = T log ag at 
which the " gas-liquid" phase transition occurs. 

Let To = min_i<^<o ^4(7). Thus, for every Tq < T < Tc, we obtain a value of the reduced 
activity of the liquid = ai/ao {ai is the activity of the liquid) that corresponds to the van 
der Waals normalization. 

Remark 3. In thermodynamics, the critical values T^, Pc, and pc are evaluated experimentally 
for almost all gases, and therefore the critical number of degrees of freedom can be set 
in advance. According to numerical calculations for a real gas, the parameter A = 1/A 
(1.6 < A < 3, Tr > 1/3) determining the point 70 ensures that the binodal passes through 
the triple point (see Sec. 3.4). The triple point can be determined experimentally with a 
sufficient accuracy. 

2.4 Holes in the Bose condensate as observable quantities. 
The maximum density of holes 

The molecules of an ideal gas can be thought of as tiny balls. Let us imagine holes, excitons 
in glass, also as balls which are empty, without the substance of a molecule. Obviously, if one 
mixes these balls in a glass in a chaotic way, then the chaos in the glass becomes increased. 
This means that the entropy increases in the presence of holes. Therefore, to achieve the 
maximum of the entropy, we must also additionally mix holes into this glass. 
In our conception, holes occur for 7 < 0. 

In the ideal gas model, we ignore the attraction, and this means that, when "stretching" 
the liquid, which results in holes, the liquid does not resist (as the sand, which is incom- 
pressible under the compression and does not resist under "tension;" cf. the appendix to the 
book [16]). 

Once there is no attraction, there is no negative pressure 'under the "tension", i.e., there 
is no formation of holes. If 7 < 0, then the plane {Z, P) is positive again, and therefore it is 
covered by the other sheet. It can readily be seen that the lines entering the point Z = 0, 
P = (i.e., to the point 7 = 0) are reflected on this second sheet back, along the same 
line. This means that it is geometrically convenient to arrange the reflection of vectors on 
the second sheet by using the matrix — /, where I stands for the two-dimensional identity 
matrix, i.e., to flip (carry out the mirror reflection for) the sheet 7 < to the negative 
quadrant. 

Note that this procedure is compliance with the concepts of Dirac hole theory, just in 
the opposite direction, namely, to a hole we assign a negative pressure, i.e., a negative 
energy. Now the straight lines can be continued through the origin to the negative quadrant, 
although the pressure really does not change its sign. This is only a convenient geometric 
"unif ormizat ion . " 

Note also that, due to absence of attraction, an ideal liquid is completely plastic; namely, 
it does not try to return to the original state (the state before stretching). In this sense, the 
Bose condensate for 7 < 0, which leads to this "kind" of liquid, can also be treated more 
visually as a glass or an amorphous solicjfl. This makes it possible to interpret the state of 
the liquid for 7 < more intuitively. 

^Physicists know that glass is a hquid and an amorphous metal is a glass. Hence, an amorphous metal is 
a liquid. It is probable that the reader will interpret excitons (holes in amorphous metals and voids in glass) 
in a simpler way than holes in liquids because the notion of holes in crystal metals is rather customary. 



Remark 3. The author has come to the revision of the thermodynamics when study- 
ing economics in which money is the very particles, according the correspondence principle 
derived by Irving Fisher. Fisher himself did not referred to his observation as the correspon- 
dence principle. However, since he was a disciple of Gibbs, there is a clear reason for the 
fact that the relation of the basic law of economics 

PQ = Mv, (37) 

where Q stands for the amount of goods, M for the number of money, v for the turnover 
rate, and P for the price of goods, is obviously related to the correspondence of economical 
and thermodynamical quantities, namely, the volume V corresponds to the amount of goods 
Q, the number of money M to the number of particles N, the rate v to the temperature T. 
The price of goods P is related to pressure to a lesser extent; however, it is denoted by the 
same symbol. 

In this correspondence principle, it is natural to correspond holes to debts and acquitting 
to annihilation. 

As mentioned above, the locus on which the chemical potential is zero gives the points 
of maximum entropy. We refer to these points as the "new spinodal." In economics, this 
new spinodal means a kind of limit for debts [151 HT] . 

Thus, according to the relations thus obtained, we face a double covering of the plane 
{Z, P} for 7 > and — 1 < 7 < 0. The meaning of the second sheet is that, for — 1 < 7 < 0, 
the chaotic state of liquid (as a phenomenon associated with the Bose condensate) increases 
when the number of holes of the type of Frenkel excitons increases, and the holes are placed 
in the liquid, which is fluctuationally concentrated on a rather slow-moving domain {from 
the point of view of the device discussed abov^, in the form of chaotic nanoholes, then the 
structure of the liquid becomes chaotically stretched. 

Here the holes-excitons cannot be indexed by our device, as well as the particles, and 
we can speak only of the density of holes. As was already said above, it is more convenient 
to place the second sheet under consideration in the quadrant [— Z, —P], by continuing the 
straight lines (19) through the singular point of Z = 0, P = to the negative quadrant. In 
other words, to make a reflection with the help of the matrix —E, where E stands for the 
identity matrix. 

Thus, it becomes convenient to speak of " negative pressure", although we neglect the 
attraction of particles, and hence there can be no negative pressure at all. As a rule, the 
pressure, as well as the temperature, is regarded as a positive quantity. We stretch the liquid, 
and it becomes plastically frozen up in this stretched state and does not tend to shrink back. 

Let us explain from the point of view of physics why the extension to the negative square 
is natural. We compare the new ideal liquid with sand, which is incompressible under the 
"compression" and "dost not resist" under stretching, because there is no attraction between 
the grains. 

Example 3. Consider a cylindrical vessel, filled with sand, whose lid is attached to the 
piston, in the room of the space station. The increase in the vessel with the piston leads only 
to a rearrangement of sand and its transformation to a fioating " dust" in the new volume 
(see m)- 

If we take into account the gravitational attraction between the grains, then the phe- 
nomenon of pulling the piston creates a negative pressure, and thus it is natural to pass to 

^In reality, the holes can change places with each other and with holes in the surrounding gas quickly 
and imperceptibly for the eyes and for the device. 
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the negative quadrant on the (P, Z) diagram, and then to neglect the gravitational attrac- 
tion. 

Neglecting the presence of attraction here is just as "legitimate" as it is in the theory of 
vapor-liquid equilibrium, where the condition that the pressures are equal is possible only if 
we neglect the surface tension. 

This also explains a smooth transition (without a phase discontinuity of the first kind) 
of this structure into ice, namely, a frozen glass crystallizes. 



2.5 Critical exponents as observable quantities under the Wiener 
quantization and the derivation of the Maxwell rule 



Mishchenko and the author [19] considered the transition to a two-dimensional Lagrangian 
manifold in the four- dimensional phase space, where the pressure P and the temperature T 
(the intensive variables) play the role of coordinates and the extensive variables (the volume 
V and the entropy S) play the role of momenta for the Lagrangian manifold, where the 
entropy is the action generating the Lagrangian structure. 

Seemingly, there is no global canonical transformation leading to a change of this kind. 
This does not confuse physicists. For example, in §25 of [2J, "Equilibrium of a solid in an 
external field," it is said that "from the equation 

dE = TdS + fidN, (38) 

represented in the form 

dS=^- ^dN, (39) 



we see ..." 



However, formula (38) does not imply the expression "represented in the form" (39). 
Nevertheless, this "implies" the following conclusion: "If the field is absent and both /i and 
T are constant, then the pressure is automatically also constant." At the same time, the 
same textbook states that, at a temperature slightly below the "dew point," "when the 
radius of the drop becomes greater than the critical value, it can be seen that the pressure 
of the liquid inside the drop differs from the pressure in the saturated vapor. The external 
field is absent. Is this still thermodynamics? Other words are used; one speaks of a vapor 
instead of gas and of the process of nucleation instead of the vapor-liquid equilibrium. And 
then a patch is immediately put on the same hole, namely, an extra term is added to (39). 
The old thermodynamics has many patches of this kind. 

It turns out that this complex transformation, leading to relation (39), can be carried out, 
as we have seen, only by continuing to the domain of negative energies. After this, one can 
justify the Maxwell transition by introducing a small dissipation (viscosity). The introduc- 
tion of an infinitesimal dissipation enables one to simultaneously solve the problem of critical 
exponents, without using the scaling hypothesis, on which the method of renormalization 
group is based. Let us show this. 

In thermodynamics, the viscosity is absent. However, generally speaking, without an 
infinitesimal dissipation, an equilibrium in thermodynamics should not be attained. There- 
fore, it is natural to implement the occurrence of this infinitesimal viscosity and then pass 
to the limit as the viscosity tends to zero. 

The geometric quantization of the Lagrangian manifold (see [20], §11.4) is usually 
associated with the introduction of the constant h. The author introduced the term of Wiener 
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(or tunnel) quantization to describe the case in which the number h is purely imaginary 

[SUES]. 

Let us apply the Wiener quantization to thermodynamics. The thermodynamic potential 
G = fiN is the action S = j pdq on the two-dimensional Lagrangian manifold in the four- 
dimensional phase space qi,q2,Pi,P2, where qi and q2 are the pressure P and the temperature 
T, respectively, pi is equal to the volume V, and p2 is equal to the entropy of S taken with 
the opposite sign. All other potentials, namely, the internal energy E, the free energy F, 
and the enthalpy W are the results of projecting the Lagrangian manifold to the coordinate 
planes pi,p2, 



E = -Jqdp, g = {91,92}, P={pi,P2}, 

W = - {q2 dp2 + 91 dpi), F = {qi dpi - ga dp2). (40) 



Under the Wiener quantization, we have 

N ^ V ^ S ^ 
9/i ' dp' dT 

Consequently, the role of time t in the quantization, is played by log(— /x/T), 

/i (9 d 



G = fiN e 



Td{ix/T) d\og{- f^/ry 



Note that the tunnel quantization of the van der Waals equation (vdW) as e — )■ gives 
Maxwell's rule (see below). 

As we shall see below, the critical point and the spinodal point are focal points, and 
therefore, as e — )■ there points do not come to the "classical" picture, i.e., to the van 
der Waals model. The spinodal points, which are similar to turning points in quantum 
mechanics, can be approached by the Airy function, whereas the critical point, which is the 
point at which two turning points are generated (two Airy functions), can be approached 
by the Weber function (see I23j). It is the very Weber function which is used to express 
the creation point of the shock wave for e — )■ in the Burgers equation is expressed. If one 
passes to the limit as e — )■ outside these points, then we obtain the vdW-Maxwell model. 
However, the passage to the limit is violated at these very points. Therefore, the so-called 
Landau " classic" critical exponents [2] drastically differ from the experiment. The Weber 
function give singularities of the form whereas the Airy function gives a feature of the 

form 

Let us present a more detailed consideration of the Burgers equation. 
Consider the heat equation 

du e d^u _ , , 

where £ > is a small parameter. As is known, all linear combinations 

u = XiUi + X2U2 (42) 

of solutions Ml and U2 of equation (41) are solutions of this equation. 
Let us make the change 

u = exp{—w{x,t)/e). (43) 
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dw 


IfdwV 


e d'^w 




' 2\dx) 
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We obtain the following nonlinear equation: 

0, (44) 

which is referred to as the integrated Burgers equatioij^. Obviously, to any solution Ui of 
equation (41) we can assign a solution Wi = —elogUi of the equation (44), i = 1,2. To the 
solution (43) of equation (41) we assign a solution 

w = —e log (e E + e = j 
of the equation (34) where /x, = — elogAj, {i = 1.2). Since 

limw = mm(wi,W2), 

e-s>0 

we obtain the (min, +) algebra of the tropical mathematics 

To find solutions for t > tcr, Hopf suggested to consider the Burgers equation 

dv dv £ d'^v 

m^'d^-2d^ = ^^ -U=Po(x), (45) 

and to refer to the function pgen = lime_j.o v (Riemann waves) as a (generalized) solution of 
the equation 

dp dp 

di^^d^^^' Pk=o = Po(a;). (46) 

The solution v of the Burgers equation can be expressed in terms of the logarithmic 
derivative 

8 

V = —e—\ogu (47) 
dx 

of the solution u of the heat equation 

du ed'^u ( 1 r "1 

m = 2d^^ = 1 - ^ j_Joix)dxj. (48) 

Thus, the original problem reduces to the study of the logarithmic limit of a solution of the 
heat equation. As is known, the solution of problem (48) is of the form 



u = {2netr'/^ exp { " ((^^ " 0' + 2t PoiO j 2t/i| di. 



(49) 



The asymptotics of the integral (49) can be calculated by the Laplace method. For t < tcr, 
we have 

u = (|J|-^/'(e(x,t),t) + 0(6)) exp ( - -Six,t)]. (50) 



Here 

.Tit) 



r-ryt) 

S{x, t) = pdx, 

J —oo 



and the integral is evaluated along a Lagrangian curve A*; r(x) is a point on A*. For t > tcr, 
there are three points ri{x), r2{x), and r^lx) on A* whose projections to the x axis are the 



^The usual Burgers equation can be derived from equation (44) by differentiating with, respect to x and 
by using the substitution v = dw/dx. 
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same; in other words, the equation Q{t,C,) = x for x in{xi,X2) has three solutions C,i{x,t), 
^2ix,t), and 



Write 



/r{x) 
pdx for X < xi, X > X2, 
-oo 

S{x,t) = min(5i, ^2, iSs), and 



Sj = I pdx, 



oo 



where J G {1,2,3} for x G [xi,X2]. 

These arguments enable us to obtain a generalized discontinuous solution of (41) for the 
times t > tcr- It is defined by a function p = p{x,t) defining the significant areas [21] of 
the curve A*. Note that this, in particular, this implies the rule of equal areas, which is 
known in hydrodynamics for finding the front of a shock wave whose evolution is described 
by equation (36). Note that this precisely corresponds to the Maxwell rule for the vdW 
equation. 

The solution v = v{x, e) of the Burgers equation at the critical point x = p^ is evaluated 
by the formula 

As X — )■ 0, after the change -4^ = rj, we obtain 

v{e, x) — T-aj^o ■ const. (52) 

What does this mean in terms of classical theory and classical measurement, when the 
condition referred to in the book |T] as the " semiclassic condition" is satisfied (i.e., for the 
case in which we are outside the focal point)? For the Laplace transform, this means that 
we are in a domain in which the Laplace asymptotic method can be applied indeed, i.e., in 
the domain where 

1 / pa: — S(p) 



u{x) = — I e ^ dp. (53) 
If the solution of the relation 







(54, 

is nondegenerate, i.e., 7^ at the point at which ^ = x, then the reduced integral (53) 
is bounded as e — )■ 0. For this integral to have a zero of the order of e^^^, we must integrate 
it with respect to x after applying the fractional derivative D~^^'^. The value of as 
applied to 1 (the value of Z}~^/^1 is approximately equal to x^^'^. 

By pSlITT] . the correspondence between the differentiation operator and a small param- 
eter of the form D ^ l/e is preserved for the ratio —e ^"/'^^ ^ while the leading term of the 

asymptotic behavior is not cancelled in the difference between ^ ^^^^ and (iue to 

the uncertainty principle (see Remark 4). 

Remark 4. Let us repeat the calculations in pT] with regard to the fact that, on this 
class of functions, D has the properties / ipD{p dx = ^ J Dip^ dx = and J xy?^ dx = 0. 
Consider the obvious inequality 



dtp 

axip + — — 
dx 



+00 

dx > 0, (55) 
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where a is an arbitrary real constant. When evaluating this integral, we have 




We obtain 

a\Axy -a+ -4|AD|2 > 0. (57) 

For this quadratic trinomial (in a) to be positive for all values of a, it is necessary that the 
following condition be satisfied: 



4(Ax)2— |AD|2 > 1 

or 

^/(K^^]KE^ > ^. (58) 

Thus, the tunnel quantization explains both the condition /x = for photons and the condi- 
tion /i < for bosons. 

In the case of thermodynamics, the role of x is played by the pressure P, and the role of 
the momentum p is the played by the volume V. Therefore, V ~ P^/^, i.e., 

~ (^ - V,y. (59) 

This is the very jump of the critical exponent. One van similarly obtain other critical 
exponents (see [25]). For the comparison with experimental data, see the same paper. 

Unfortunately, thermodynamics does not use the concept of Lagrangian manifold which 
was introduced by the author in 1965 [26]. It is especially suitable for thermodynamics, 
in which there are pairs of intensive and extensive quantities. Intensive quantities, roughly 
speaking, are the quantities for which one cannot create the concept of "specific" quantity. 
These are the temperature T, the pressure P, and the chemical potential fi. To these intensive 
quantities, there correspond related extensive quantities, namely, the entropy S, the volume 
V, and the number of particles A^. Altogether, they form the phase space, where the role of 
coordinates is played by the intensive quantities and the role of momenta is played by the 
extensive quantities. In this case, a Lagrangian manifold is a three-dimensional submanifold 
(of the six-dimensional phase space) on which there is an action, an analog of the integral 
S = J pdq, q G M^, p G in mechanics. It is locally independent of the path. 

Usually a 4-dimensional phase space T, S; P, V is considered. This space corresponds to 
q G R^,p G R^, qi T, q2 P, Pi S, p2 ^ V: dS = pi dqi + p2 dq2, depending on the 
coordinate plane of the form qi,q2', qi,P2', Pi, 12] Pi,P2, to which the Lagrangian manifold is 
projected. 

The Lagrange property means that the number of planes cannot coincide (there are no 
planes of the form and q2,P2)- To every projection there corresponds some potential 

(gi,g2 is the thermodynamic potential, etc.). 

This is an obvious correspondence. If it were more elaborated, then, on one hand, the 
transition from the action dS = pidqi +P2dq2 to the "action" -coordinate dqi = dS/pi — 
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iPi/vi) dq2 would be not "so obvious" (see formulas (28)-(29)). On the other hand, it would 
be natural to use the semiclassical (Wiener-Feynman) quantization of action rather than 
the scaling hypothesis. 

(The quantization of the Lagrangian manifold differs from the full quantization of ther- 
modynamics [27] in the same way in which the semiclassical Bose-Sommerfeld geometric 
quantization differs from the quantization of Schrodinger, Heisenberg, and Feynman. 

The term "dequantization," which is well-known in the tropical mathematics [24], means 
the Wiener or tunnel quantization of the Lagrangian manifold, and then the passage to the 
limit as the quantization parameter (the viscosity) tends to zero. 

3 Considering the attraction. 

Dimers (pairs) as observable quantities 

3.1 Second quantization of classical mechanics and ultrasecond 
quantization of thermodynamics. Operators of creation and 
annihilation for pairs-dimers 

The second quantization is always associated with the identity of particles, and, if it is carried 
out for classical particles, then it is tacitly assumed that the particles are indistinguishable 
for the observer. Instead of an A^-dimensional problem, we arrive at the three-dimensional 
picture in which N particles are distributed. The Vlasov equation [281 129] is obtained from 
the second quantization of classical mechanics. However, the original arguments used by 
Vlasov were actually based on the assumption that the particles can be regarded as identical 
ones. 

In the classical system of Hamilton equations for particles, even if the Hamiltonian is 
invariant under any permutation of the particles, the initial conditions need not have this 
property. However, the initial conditions of the Liouville equation can be regarded as data 
satisfying the conditions of symmetry. 

Indeed, let 

p°=p°i,p°2,py (60) 

are initial conditions for the Hamiltonian system whose Hamiltonian is invariant under every 
permutation Pi,Pj and also under every permutation qi,qj. For example, let 

^(^'^) = Ei^ + EE^(i^^-^.i)- (61) 

Substituting the initial conditions into the Hamiltonian (61), we obtain the energy 

E = H{p\q'). (62) 

The energy is conserved along the trajectories of the Hamiltonian system. 

Consider further the Liouville equation corresponding to the Hamiltonian system 

|^ = {if(p,g),p}, (63) 

where {■} stands for the Poisson bracket, with the initial condition 

p(E), p{E) G 

5>fi 



where E satisfies (62). This equation describes the distribution corresponding to the Hamil- 
tonian system with the initial conditions given above. 

The Liouville equation and the initial conditions are symmetric with respect to any trans- 
position of Pi and pj and to any transposition of and qj. This symmetry is preserved for 
the solutions. According to the Gibbs distribution for the Gibbs Ensemble, every distribu- 
tion can be expressed in terms of energy. Therefore, it is symmetric with respect to any 
permutation of the particles. 

Schonberg [30] carried out a second quantization of this system in the Fock spac^. 
In [29], the Vlasov equation was obtained under the assumption that the interaction is 
small and the number of particles is large. The BBKKI chains are also symmetric with 
respect to these permutations. Hence, for any distribution in the many-body problem, such 
a symmetry follows. Thus, we arrive at the invariance with respect to the permutations of 
the particles, and thus to a "distribution of Bose-Einstein type" for the statistical physics of 
classical particles. Other mathematically rigorous arguments which lead to a "distribution 
of Bose-Einstein type," in the form of lemmas and theorems (see [31]), and hence also the 
distributions of classical particles, obey the laws of number theory. 

Although a modern macroinstrument cannot trace the motion of every particle (because 
of discreteness of the observation times in relaxationally stepwise process), which is possi- 
ble at the classical level, but it can distinguish between molecules and dimers or clusters. 
The clusters consisting of more than two molecules occur in gas much less frequently than 
dimers. Dimers are observable at all temperatures, and a macroinstrument can calculate 
their average percentage at a given temperature. This is an important new phenomenon in 
experiment, and this phenomenon was not available to the great who formulated the basic 
laws of thermodynamics. 

Dimers occur and become immediately split by monomers (single molecules), and they 
are created and annihilated in different places. They occur because there are quantum forces 
of attraction between molecules (the dipole-dipole interaction). The dimers are virtual, as 
the ideal liquid is. 

To take into account this important phenomenon (the creation and annihilation of dimers) 
mathematically, one obviously needs to make the "second" quantization and introduce the 
creation and annihilation operators for dimers, i.e., for pairs. The author referred to this 
"second" quantization as the "ultrasecond" quantization due to the introduction of creation 
and annihilation operators for pairs. In the special case of the Bardeen model, this quanti- 
zation was introduced in essence [32], [21]; however, since this model is, roughly speaking, 
exactly solvable, these operators turned out to be hidden in the model in a sense. 

In general, the ultrasecond quantization and the asymptotic behaviors associated with 
it are rather cumbersome and lead to quantum equations involving the Planck constant h. 
The passage to the limit as ^ — )■ 0, and then the passage to the limit as the viscosity tends to 
zero, are cumbersome, and we present here only a part of this passage, which is related, as 
in the previous section, with the introduction of an infinitesimal viscosity into the classical 
scattering problem (for ^ = ). This means that we introduce the viscosity and, after 
manipulations, pass to the limit as the viscosity tends to zero. This procedure will enable us 
to find the Boyle temperature Tg, and then also the Boyle density pb, i.e., the so-called Zeno 
line, which is present in the van der Waals model and which was first noticed by Bachinski 
in experiments with pure gases. 

We shall further obtain the so-called law of corresponding states. 

^This space exists only under the assumption that the particles are indistinguishable (the "commutativ- 
ity," or the invariance). 
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3.2 Boyle temperature as the temperature above which the dimers 
are not observable in the Boltzmann— Maxwell ideal gas 



The attraction between the particles occurs in the quantum mechanical consideration of the 
dipole-dipole interaction. In the standard semiclassical limit, if the distance between neutral 
molecules is fixed (does not depend on the parameter H, i.e., on a dimensionless parameter 
proportional to K), then, as /I — > 0, the attraction disappears. In this sense, the use of an 
attraction potential in molecular dynamics using the classical Newton equations for many 
particles is at least baseless. 

However, if, along with h, the problem involves other small and large parameters, then 
the attraction potential can be kept for some relationships among these parameters under 
the passage to the limit as /I — )■ 0. 

Since the scattering problem has another parameter tending to infinity, for example, the 
time of scattering is considered in the interval from — oo to +oo, it can happen that, ash ^ 
and t — > oo simultaneously (provided that there is a dependence between these parameters) , 
an attractive potential of the order of is kept (as h ^ 0), where r stands for the distance 
between the particles. 

As we see below, to obtain a "rough" thermodynamics leading to the law of corresponding 
states, it is sufficient to determine the values of the Boyle temperature Tb and the Boyle 
density pb only for mercury. As is known, the mercury isotherms are very close to the Van 
der Waals model (see Fig. 14), and hence the Lennard-Jones potential model must provide 
a good description of the following two important facts: 1) attraction existence; 2) collision 
of molecules due to a rapidly increasing repulsive potential. 

We present only a typical example of studying the relation between the actual gas and 
the interaction potential, which corresponds to the case of a small intermolecular distance 
such that quantum effects must inevitably arise. But the natural choice of the potential Q 
for an actual gas and and the fact that the Zeno lines are taken into account give us mercury 
isotherms in Sec. 3.5. 

As an example, we consider the Lennard-Jones potential, noting that, in our fundamental 
manipulations, the repulsive part of the potential does not play any role. 

The only essential quantity is the so-called effective radius a, because it determines a 
one-dimensional elementary length. 

As is known, in the radially symmetric case, 

mv'^ ^ , , ^ , 
+ ^ + $ r = E. 64 

In the original scattered particles, we prescribe an energy E and an impact parameter 
B. The momentum M, as well as the energy E, is preserved. We also know that 

= B^E. (65) 

Expressing the energy E, we obtain for the attraction 

(mt>^)/2 + $(r) 
^ = 1 - Byr^ ^^^^ 

in the domain r < B. 

In the scattering problem, for the interaction potential, one considers the Lennard-Jones 
potential 

^(^''0=4.(p^^-p^^), r = r'-r", (67) 
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where e stands for the energy at the well depth, a for the effective radius, and ||r' — r"|| 
for the distance between two particles with radius vectors r' and r" . In the two-particle 
problem, in the absence of external potential, the problem is reduced to a one-dimensional 
radially symmetric problem. 

In problem (66), for different values of B, there are other barriers and wells. At the 
stationary points -Emin and -Emax, the velocity vanishes, and thus these values can be evaluated 
by using the potential term only. 

We speak now of a pair of particles with the mass center which is caught by the trap 
(rather than of a single particle). Therefore, the difference -Emax — -Emm is the very energy 
which is needed to knock out the pair (dimer) from the trap. 

After formation of a dimer at T < T^, one should consider the collision of a dimer with 
a monomer according to the same scheme, assuming that the pair-dimer (a yU-particle) has 
already been formed. Further, one considers the scattering problem already for the dimer 
at -Emin = T and a particle of the same mass m. The successive consideration of such 
a multistep procedure leads to formation of a three-dimensional cluster, and we obtain a 
temperature significantly lower than Tc. It is important that the first step gives an upper 
bound, and this restriction is natural. 

In an experiment, the percentage of dimers in gas can be calculated. It can be seen how 
dimers are created and how they are annihilated (broken by monomers). After this, the 
average number of these events is calculated. The higher is the temperature, the higher is 
the average energy of the monomers, and the smaller is the number of dimers. 

The main point is that, under this approach, there are only two values, -Emax and E'min, 
which are kept in the skeleton of the scattering problem (cf. the skeleton of the amoeba in the 
tropical mathematics P^]). For E^max = E^min, the well disappears. For the attractive part 
of the Lennard- Jones potential, this energy is equal to O.Se. With regard to the standard 
Clausius considerationqj, we can see that the average energy of the particles is equal to ^e. 
The average energy is the temperature, and T = y|. Above this temperature, there is no 
well. In thermodynamics, for physical reasons, this is the so-called Boyle temperature Tb- In 
our framework, the Boyle temperature is defined as the temperature above which the dimers 
are practically absent. This is a new approach. According to this conception, the Boyle 
temperature for argon (Ar) is Tg = 382 K and for krypton (Kr) it is Tg = 547 while the 
tables of the experimental work |33] give Tg = 392 K for Ar and Tg = 538 K for Kr. The 
discrepancy between the theoretical and experimental values is of the order of 2-3%. 

The critical temperature -Emax must correspond to the deepest well, i.e., to the maximum 
value of the difference E^max — E'min for all impact parameters B. This difference determines 
the drop of the energy of a dimer after this dimer was captured by the "trap," and thus 
determines the energy which a monomer must have to knock out the dimer from the well 
(i.e., for the dimer to collapse). 

The height of the barrier "protects" the created pair whose reduced mass was captured 
by the trap of "dimers" and "clusters" from "shocks" of monomers. As the temperature 
decreases, T < Tc, the height of the barrier reduces, and, to survive, the clusters must create 
their own barrier in the form of a microanalog of the surface film. Thus, a "domain" must 

^Following Clausius, experts in molecular physics usually argue by proceeding from the symmetry of the 
motion of a molecule in all six directions. In the scattering problem, we use the principle of symmetry in 
all directions, which is standard in molecular physics. The fraction of all particles that moves head-on is 
1/12. There are three such directions; hence, one quarter of all molecules collide. The arguments concerning 
symmetry that were used by Clausius to evaluate the free path length (and are repeated here by the author) 
are quite approximate. However, these arguments do not modify the values of the ratios of the form Tg/Tc. 
This very ratio is of interest for us. 
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occur, a three-dimensional cluster (the so-called elementary cluster) which has at least one 
particle which is protected by other particles. 

This is a new definition of the critical temperature Tc as the temperature below which 
clusters are formed from dimers. Calculations give E ^ax = 0.286| at the point maxB(_E'max — 
-E'min)- The impact parameter at this point is equal to B = 2.436a. 

3.3 Macroinstruments and microinstruments in dimension theory, 
determinin the maximum density pB and the Zeno Hne which 
borders the domain of dimers as the density is modified 

Let us now obtain analytical formulas for the Zeno line in dependence on the potential. 

E/e 
1.0 

0.8 

0.6 

0.4 



"■"l.O 1.5 2.0 2.5 3.0 3.5 r/a 

Figure 6: Wells and barriers in the scattering of two particles with the Lennard- Jones 
interaction potential at different impact parameters B. 



T 

Figure 7: A trap for a fictitious particle-dimer, in the center of mass system (CMS). Here 
r stands for the radius vector of the dimer; it is marked on the abscissa axis. A "particle" 
falls from the left, from the point r = B, where B stands for the impact parameter. After 
refiection from the "wall", i.e., from the potential, in the presence of viscosity, the particles 
meet the barrier and then after rather many reflections from the barrier and the "wall" drop 
on the well bottom in a greater time interval. 

Let us use considerations of dimensional analysis for the scattering problem and the 
definition of the one-particle (the so-called thermal) attraction potential. The scattering 
problem is considered for the Lennard-Jones potential, and therefore there is an additional 

an 





parameter of length in the problem, the parameter a, which is the effective radius. The 
attraction potential occurs in the quantum theory of the dipole-dipole interaction. If we 
fix the distance between the particles and assume that the semiclassical parameter h tends 
to zero, then the attractive potential vanishes. This means that, as /i — )• 0, the distance 
between the particles decreases. It follows that the attractive potential acts between "nearest 
neighbors" only. 

Therefore, it is natural to use the expansion of an attraction single-particle potential in 
powers of the radius r up to 0{r^/V) only. 

The "dressed" or thermal potential \I'(r^) is attracting. As is well known, it was derived 
from the fact that the correlation sphere for the A^-particle for the Gibbs distribution is 
finite [31]. 

One of the most interesting points of independence of a macroinstrument of a microin- 
strument manifests itself when applying the dimension theory P3]. A macroinstrument 
determines the volume V. According to the independence of the thermodynamic quantities 
on the shape of the volume V, the volume V ensures us that we have the dimension in 
the dimension theory is three; however, this volume does not give us any one- dimensional 
measure, which is the typical length in the thermodynamic process. 

A microinstrument determines the effective radius of the molecule and the mean free 
path. However, in dimension theory, we cannot measure the typical length of the macrother- 
modynamics by using the radius of the molecule or the mean free path, even if the volumeis 
a Torricelli tube and its typical thickness is small as compared with the case in which the 
vessel is a ball. These considerations show that the only possible dimensionless combination 
for the argument of \l/(r^) is 

•qjl — 
V V 

Since <^ V, it follows that the expansion is 

The constant Ci gives no contribution to the scattering problem, and the thermal "single- 
particle" potential turns out to be proportional to the density. On the plane {T,p}, the 
maximum of the binodal (according to Fig. 8 for T = T^) is equal with respect to p to the 
very value pc, which enables us to find the proportionality coefficient. It turns out to be 
equal to one. 

The situation in which the thickness of the tube is "commensurable" to the radius of the 
moleculj^ leads to quite different effects: to the superfluidity of water in a nanotube and to 
the freezing of water at T° = (see f36]). 

Expanding 

2 2,2 (ri-ra)^ (ri + ra)^ 
r =ri+r2 = + , (69) 

we can make the separation of variables in the two-particle problem to a scattering problem 
for a pair of particles and to the problem of their joint motion for ri -|- r2, just as it was 
done in [37]. In this case, in the scattering problem, a quadratic attraction potential (an 
upturned parabola) — pr^ is added to the Lennard- Jones interaction potential, and p = y 
for the isochoric process. 

In the scattering problem thus obtained, there are two points of rest, namely, the stable 
one, -Emiiu and the unstable point -Emax- Their ratio is a dimensionless quantity. As follows 



"'^'^When the scales become "commensurable" in this sense, another thermodynamics arises [371 [ 
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Figure 8: The T-p-diagram for gases corresponding to simple liquids. Here Tr = T/Tc and 
Pr = p/Pc- The Zeno line (the straight line Z = = 1.0) in the phase diagram. For the 
states Z > 1.0 (the hard fluids), the repulsive forces dominate; for the states Z < 1.0 (soft 
fluids), the attraction forces dominate. 

from the previous section, Z = j^, where P stands for the pressure, N for the number of 
particles, T for the temperature, and, due to the fact that a stable stationary point has the 
meaning of temperature, it follows that the ratio 



PV TP 

Z=^ = ^ (70) 
NT E 

enables us to wrote the curves Z = const in the graph T, p = y. 

The curve at Z = 1 is called the Zeno line (or the Bachinski parabola), and the locus of 
the beginnings of the curves ^max (for C2 7^ and B — )■ 00) is referred to as the binodal. 

Thus, in our view, the Zeno line determines the temperatures for which the dimers become 
practically nonexistent for a given density. 

Denote by pc the value of p at the maximal point of the binodal and denote the endpoint 
of this curve on the p axis by p^. In [38]. this point was referred to as the hypothetical 
point pb (the Boyle point). 

Calculating the value of Z^., we obtain Zc = 0.296, which coincides with the values of Zc 
for the noble gases up to thousandths. The ratio pc/ Pb also coincides with the values of this 
quantity for the noble gases. 

Table 1 shows the data corresponding to the resulting diagram (for B = 100 in "molec- 
ular" values), and note the discrepancy between the basic dimensionless relations obtained 
by the data of molecular dynamics and the theoretical relations obtained by physicists from 
the chain BBKKI and the A^-partial Gibbs distribution. 

Table 1 





PcIpb 


TJTb 


0.29 


0.273 


0.36 


0.308 


0.285 


0.38 


0.375 


0.333 


0.296 



On the top line of the table, the theoretical values for Zc, Pc/ Pb-, and Tc/Tb, obtained 
using the above theory are presented. The second line contains the values of the same quan- 
tities evaluated according to the latest data of molecular dynamics and results of theoretical 
physicists for the Lcnnard-Jones potential. The third line gives the values obtained from 
the van der Waals equation, which is empirical. 

The value of can be computed in the experiment very accurately, and it is equal to 
0.29 for noble gases, nitrogen, oxygen, and propane. The value of pj Pb (the ratio of the 
critical value of p critical to ps, i.e., to the entire length of the segment with respect to p on 
which the Zeno-line "cuts" the abscissa axis away) evaluated in the above theory coincides 
with the corresponding values for water, argon, xenon, krypton, ethylene, and a number of 
other gases. 

Let us present detailed calculations to find the Zeno-line. 
Consider the potential 

-ar"^ + r^U (r) 

Its first derivative is equal to 



_ r{2B^U{r) + r{2ar{-2B^ + r^) + {B^ - r^)U'{r))) 



and the second derivative is 



1 



E"(R) = -—±—-(2{B^ + Wh^)U{r) + r{-2ar{QB^ - 35V + r^) 

+ 4(S^ - B\^)U'{r) + r^{B^ - r^fU"{r))). (73) 



We obtain a solution of the equation in the form 



^' -8U{r) + 2rU'{r) + 2r^U"{r)' ^ ' 

Substituting the value B{r) into (66), we find E{a.), the Zeno-line, i.e., a segment, which 
is straight up to 3%, T/Tb -\- p/ Pb — 1, where pb stands for the maximal density as T — > 0. 

3.4 Limit stretching of a liquid. 

The maximum density of holes 

We have Z — P/ pT, where p = (1 — T/Tb) is an isotherm-isochore. Therefore, 

Pr 

PBTril-T/TsY ^^^^ 
Let us express Z in terms of 7 for Z < and 7 < and for /i ~ o(l/ logN), 

N = A('y)T P = — ^ ' 




Figure 8A: Curve 1 — Zmin- 

where Tr = T/Tc and Pr = P/Pc- The value of pB is defined in this very normahzation. 
Therefore, 

T,^+X(2 + 7) _ TX(2 + 7) 
C(2 + 7c)(A(7)T2 C(2 + 7c)A(7)' 

On the other hand, N = p = ps {1 — T/Tb) = A{'y)T. Consequently, Pb/T — Pb/Tb = ^(7) 
and 

/ I I \\'y\ 



(76) 



TJ \pb Tb 



Substituting this into (76), we obtain 



Z 



1 . O'^l 



C(7 + 2) Vpj; Tb 

C(7c + 2) A(7) 



and ^(7) — )■ cxD as 7 — 7- —1. 

Therefore, as 7 — )■ —1, the value l/Tg is negligible. Consequently, 



C(7 + 2) 
p^^iC(7c + 2)A(7)i-H 



^l7^-i - ^ iil^J . (77) 



Here 7c corresponds to Zc, which is the minimal value of Z on the critical isotherm for P = 1 
(see Fig, 9). 
Since 

A(7,A) = (A>-M7))'/'+^ 

where 7 < 0, it follows that 

A(7,A)i+^ = A>-M7), 

where 

c(7) = 



t"'dt{- ^ 



J e*-l 
Hence, 

C(2 + 7) 



Z 



p^^'C(2 + 7c)c(7)A>+i 

The expression 0(7) tends as 7 — — 1 to (1/2) loge:, where e stands for the lower limit of 
the integral expression for A{'jY~^^^ at 7 = — 1. 



Similarly, the expression C(7 + 2) at 7 = — 1 is equal to logs. 
Hence 

Z\ 2 
as 7 — > and Z — )■ 0. Therefore, 

Z < ^max < const 

for all values of 7. 
Hence, by (75), 

Pr = ZpBTr (^1 - < Z^J^. (78) 

Moreover, it is clear that — > as — )■ and 7 — > 0. 

The value of A determines the minimum of Pr and the maximal density of holes. 

To determine the gas-liquid transition, as in Sec. 2.3, with the correction to the Zeno line 
taken into account, we derive equations of the form (25)-(26), by normalized the activity for 
the critical isotherm. 

z 1 
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Figure 9: The picture of the law of corresponding states for molecules of different gases. 
Thin lines show isotherms for methane. Different symbols on isotherms correspond to argon, 
carbon dioxide, water, etc. The fact that isotherms of different gases are close to one another 
illustrates the empirical law of corresponding states. The theoretical isotherm (the solid line) 
does not fully coincide with the experimental one. This is an effect of the same type as the 
jump of the critical exponents. The viscosity (the Wiener quantization) smooths the sharp 
angle of the hmit isotherm (as £ — > 0). 



3.5 A coarser measuring instrument 

and the law of corresponding states 

Obviously, the rougher is the device, the greater is the dissipation parameter, and therefore 
the more important is the uncertainty principle. If our device does not distinguish some 
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Figure 10: Spinodal curves in the domain P < 0. The in dined hne is the continuation of 
the theoretical critical isotherm shown in Fig. 9 to the domain of negative values of P and 
Z. The initial point of the curve A = 1/3 (A = 1/A) is at a distance from the point P = 0, 
and the starting point A = 0.5086 (A = 1/A = 1.97) of the curve coincides with the point of 
intersection of the dashed line and the solid line. 

molecules from one another, then this is a rough instrument. It is not necessary that the 
device itself be so rough. It is sufficient to say that the experimenter does not want to 
distinguish molecules and computes the density roughly, counting all the molecules for which 
he wants to construct a "rough" thermodynamics. Hence, when speaking more precisely, 
from the point of view of mathematics, the rougher thermodynamics we want to receive, the 
more rough will be the answer due to the Wiener uncertainty principle. 

Our rough instrument does not distinguish particles not only in mixtures. Measuring 
different particles, the observer thinks that these are the same molecules, and only the 
measurements are somewhat scattered. This is even a more rough thermodynamics. It is 
referred to the fluid part, when the gas-hquid is not distinguished, clusters occur, etc. The 
critical points are on the boundary of the fluid domain, and we approach them from the side 
of the fluid domain for T > Tc. 

When considering above the phase transition, we have equated the chemical potentials 
of the liquid and the gas phases on the critical isotherm, assuming that there is no phase 
transition there. We considered the case in which the viscosity vanishes, and we obtained 
a phase transition which was not smeared. Now, when considering the problem with a 
coarser device, we no longer have any right to assume that the dissipation tends to zero. 
The uncertainty principle gives us a fairly large smearing of the phase transition; however, 
it happens on a "rough" critical isotherm which is measured by a rough instrument. The 
latter cannot distinguish now not only particles of a single pure gas but does not distinguish 
particles of different gases either. 

We have compared the Wiener quantization of thermodjTiamics with the dissipation 
resulting in a shock wave. However, the dissipation may be different in different substances, 
while we are interested now in the Hugoniot conditions for the entire mixture, and we do 
not want to distinguish stratiflcationally occurring internal shock waves within a common 
shock wave. This is particularly evident when the dissipative parameters of viscosity and 
thermal conductivity are different (see [IQ], § 95), the viscosity u is small, and the thermal 
conductivity is relatively large, x ^ z/. If the processes inside the shock wave are not of 
interest for us, then we pose the Hugoniot conditions on a shock wave spread with respect 
to heat conductivity. 

On the other hand, if the values of viscosity in a mixture of different gases are different but 
not dramatically different from one another, then the width of a shock wave thus obtained 




is defined by the average viscosity. 

Since, as a rule, the Bose-Einstein-type distribution is considered in the three-dimensional 
case, it would seem to be natural to add the integration over the coordinates to the integration 
of the momenta. Then the three-dimensional volume V would occur as if this is a natural 
way. It would seem that there can be a generalization if the dimension of volume is changed 
when changing the dimension of the momenta. However, in the manual [2] , when considering 
a photon gas, the integration of the Hamiltonian of the oscillator is carried out both over 
the momenta and over the coordinates, and, nevertheless, the three-dimensional volume V 
is taken as the multiplier for the distribution. 

Note first that number theory gives, for dimension two, a distribution without the volume 
V (see Example 1), as well as the initial relations in [2] both for the Boltzmann-Maxwell 
ideal gas and for the Bose-Einstein ideal gas (see (l)-(4)). Further, in [2], both the distri- 
bution for bosons and the distribution for photons are multiplied by the three-dimensional 
volume V. Certainly, the main distribution is the distribution without the volume V, and 
its multiplication by V is caused by the correspondence with thermodynamics in which the 
pair 'the volume V - the pressure P' is the most important tool. 

Therefore, the most natural generalization to the nonideal distribution is the multipli- 
cation of a fractional Bose distribution by a function of V of the form Vip{V/VQ), where Vq 
stands for some reference volume and the function (p{x) is smooth. 

The introduction of this multiplier does not change the distribution caused by number 
theory [H] in which the variable V is eliminated by the change of the variable N/V = p. 
However, if we consider V(p{V/Vq) as a multiplier, where Vq is some typical volume, then 
the relation for the three-dimensional Lagrangian manifold in the six-dimensional phase 
space {P, V; T, S; fi, N} is preserved, whereas the variables V and N do not convert here into 
a single variable p = N/V. Therefore, the volume and the number of particles are changed 
on the isochore p = const in general. This modification of the f2-potential does not change 
the specific entropy, which is also of importance. 

In a mixture of gases, we are to choose a reference gas in which the difference between the 
vapor and the liquid is the lowest possible, for example, from the point of view of the number 
of dimers. This gas is the mercury vapor {Z^ = 0.4). Let us carry out a normalization of 
activity (26) for an isotherm of this gas at Z = 0.4 and assume that there is no phase 
transition on the critical isotherm of mercury. In accordance with what was said above, the 
rough device cannot distinguish among molecules of / distinct gases. Let us calculate the 
average degree of freedom for these molecules by taking the arithmetic mean of the values 
of the entropy of I pure gases (see (13)) on the basis of experimental data for Z*, i = 1, 2, 
. . . , /, pure gases (see (13)), 

_|_ o\ Cilaverage + 2) _ i 0\ — }_ \^(^ _|_ O^l ^^i^LZL^l 

\ I average ' ^) t./ , i \ ^averaqeyfaveraqe ~r -^J , / ,\fi ~r ■^J , i \ • 

CKlaverage + 1) / ^ C(7i + 1) 

The highest value Z^. = 0.4 is given by mercury (Hg), and therefore the average number of 
degrees of freedom of this family of molecules is certainly less, and therefore Zaverage < -^c 
(mercury), 7 < •yng- The critical pressure is greater than that for mercury, CilHg + 2) < 
C{laverage + 2).. Therefore, the value we have chosen for mercury, = 1, is less than 
Paverage = (^average + 2)/({'yHg + 2). Thus, the valuc P,. = 1 for Zng bclougs to the domaiu 
of the phase transition "gas- liquid" for Zaverage- This implies that, for P,. = 1, the phase 
transition to liquid occurs at Zaverage- 

This phase transition to liquid can be seen in Fig. 9, in the form a vertical bounded by 
a black sloping line depicting the liquicf"]. 



^^One can rigorously prove the existence of phase transition only for the transition of a new ideal gas into 



Taking into account the Zeno line influences the form of the Q-potential as follows: 

poo ,7+1 7, 

(79) 

where y — exp{ii/T) is the activity and /i stands for the chemical potential. 

Let us write out the differential equations for (fi{x) with regard to the relations on the 
Zeno line, 

T, = Ts (1 - pJpb) , P. = p,Tb (1 - p,/pb) , (80) 

where the subscript z means that the corresponding values are taken on the Zeno line, i.e., 
for Z = 1. 

Let us construct the relation Z = 1 on the Zeno line. This relation is of the form 

^ ^ dn/dv ^ ifiyjvo) + iyjv^W iyjVo) Li^+2(i/.) 

It follows from (75) on the Zeno line that 

Tl+^[v{VjVo) + {VjVoW{VjVo)] U,+2{y.) = (iV/K) Tb (1 - N/{V,pb)) . (82) 

Assume that the conditions N/V = const and N = const hold on the isochore and on the 
Zeno line defined by relation (80). It follows from (81) and (82) that 

(1 - N/{V,pB)y-'' viVjVo) Li^+iiVz) = N/{V,pb). (83) 
Here N/V^ ~ const and = const, and hence this equation contains unknowns y^, V^, and 

After finding the value yz as a function of Vz and Lp{Vz/Vo), wc substitute it into for- 
mula (76), thus transforming equation (76) into a differential equation for the function ip{x) 
depending on the constants pb and Tb- The equation for (p{V) enables one to find the point 
yz{p) on the Zeno line, and, after this, the isochore is changed only at the expense of the 
modification of activity a — e"'^^'^ from yz{p) to a = 1 for the polylogarithm Lij+i{a). The 
function '^{V) is not reflected in the structure of the isochore, and it is reflected in scaling 
only (see Fig. 11, (a)). 

The most important problem in the theory of differential equations, the existence problem 
for a solution, remains open. Physicists pay less attention to this problem than mathemati- 
cians. 

What are conditions for the existence of a solution of equations (81)-(82)? 
Introduce the notation Tg = Tstandard = T/Tm and Pg = Pstandard = P/Pm, where and 
Pm are defined by the formulas 

P 7^7+2 T/ T/ T/ 

m -'- m 1 ^0 'max; 

and = ^max and Pm = Pmm are defined below. 

Since Nc does not depend on T, it follows that V and N are constant along the isochore 

Let us write out the above relations at the point = 1, = 1: 

[v{Vz/Vo) + {Vz/Vo)v'{Vz/Vo)]aj + 2) = l, K„v^(l)C(7 + l) = Ar„ = K.ax, (84) 

a new ideal liquid without taking into account the Zeno line, which is unknown for 'Yaverage- Therefore, a 
rigorously proven transition from Z = 0.4, P = 1 to a liquid is obtained a bit higher than at the point 
Z = 0.12, P = 1 in Fig. 10. 



Hence, 

_ Kn C(7 + 1) Li^+2(1) _ I , . 

iV,C(7 + 2) ■Li,+i(l) iV, p^- ^ ' 

Since we construct isochores V = const on the plane {Z,P}, it follows that Vm — Vz- 
Eliminating (V"^99(V"^/Vo)' from (84) by using (83), we obtain 

Since Kn = V^, we see that 

99(1) = iV,(C(7 + 1)K.)"\ ^(1) + ^'(1) = (C(7 + 2))-\ (87) 
and it follows from (81) on the Zeno line that the following equation holds: 

Li^+2(|/^) K« C(7 + 1) 



Li,+i(y,) iV, C(7 + 2) ^' ^^^^ 

Eliminating y^, we find a relation for and 7. The maximum value of Z„ depends on 
the values of ps and Tb only. For mercury, this maximum is obtained for 7^ = 7inin = 0.1, 
and = 0.4, which coincides with the value of Zc for mercury. This coincidence, which 
depends on pB and Tb, holds for mercury Hg only (among all the elements of the periodic 
table), which confirms the correct choice of the Q-potential (79). 

The family of isochores, according to system (82)-(83) with the above initial condi- 
tion (87), is shown in Fig. 11 (b). 

The first relation for the hmit isochore Vm/Nc = Zm, for 7 = 7^, is of the form 

7 7 C(7+ 1) ^h+^jy) 1 / / ( \ (Qa\ 

Z^Z^- ^ • T ■ . X , l<y<2/47)- 89 
C(7 + 2) Li^+i(2/) 

Since P, = Li^+2(l/)/C (7 + 2) and = T^+^ U^+i{y)Vm, it follows that the other 
relation is 

Equations (89) and (90) give an almost straight segment of the isochore. 

Starting from Z < 0.4 (for example, for a van der Waals gas), the phase transition to a 
liquid occurs for indistinguishable particles of the law of corresponding states. This gives a 
broad area (at the expense of the uncertainty principle of a "rough device" ) around the line 
segment P = 1, 7; = 1, Z < 0.4. 

Remark 5. Since the rightmost isochore (which is not shown in Fig. 11) is a segment 
of a straight line, it follows that all isochores of high density must also be line segments. 
They pass through the point p > pm on the Zeno line and the point = 1/p on the line 
Pg = 1. We thus obtain (when including the isochores shown in Fig. 11) a complete family 
of isochores for Z < 1, P > 1, related to the law of corresponding states. To any point of 
an isochore in the plane {Z, P} there corresponds a point of temperature, and we construct 
isotherms which, up to the Wiener uncertainty principle, correspond to the experimental law 
of corresponding states (Fig. 12). 

In Fig. 13, the graph of an experimental isotherm for mercury is presented; this graph 
was kindly evaluated by Professor V. S. Vorob'ev, according to the most recent data, at the 
instance of me. Note that the passage gas-liquid happens at T = 1473 K along a slanting line 
rather than a vertical one, which is related to nonzero viscosity and the Wiener uncertainty 
principle. This effect is of the same nature as the jump of critical exponents and "thickness 
of the layer" of a shock wave. 
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Figure 12: The dotted line shows the Zeno hne Z = \. The bold line is the critical isotherm 
of a real gas (mercury) which is calculated theoretically, and the fine lines are isochores of 
mercury for T < Tc. 
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Figure 13: The thin sohd hnes represent the experimental isotherms for methane (see Fig. 9). 
The dotted hnes approximating the experimental curves are based on theoretical data. The 
dashed curves show the experimental isochores. 




Figure 14: The Hougen- Watson diagram for mercury. Experimental data (from the sources 
W. Gotzlaff, G. Schonherr, F. Hensel, Z. Phys. Chem. Neue Fol. 156 219 (1988) and 
W. Gotzlaff., Ph. D. Thesis, University of Marburg, 1988) are equipped with symbols. The 
thick lines correspond to the van der Waals equation for the related temperatures. 
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